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Abstract 

The coupling constant of — > KK decay is calculated in light cone QCD sum 
rules. The result obtained for g^KK = (4.9 db 0.8) is in a good agreement with the 
existing experimental result. 
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1 Introduction 



Light scalar mesons constitute a remarkable exception of the quark model systematization 
of mesons and their nature still need to be unambiguously established [1]. 

Particularly, the nature /o(980) meson is under debate. According to the naive qq 
picture and strong coupling with kaons, /o(980) can be interpreted as a pure ss state [2]- 
[4]. However, this interpretation does not explain mass degeneracy between / (980) and 
isovector ao(980), which is interpreted as a (uu — dd)/\/2 state. It is also interpreted as 
a four quark qqss [5] bound state of hadrons [6] -[8] and as a result of a process known as 
hadronic dressing [2, 9]. 

For understanding the content of the /o meson several alternatives have been suggested: 
For example, analysis of — > f j decay [5]— [10] and investigation of the ratio T(a — > 
/o7)/r(0 — > /q7) [7, 8] are believed to be the most promising ones for this purpose. 

The — > /q7 decay is a very efficient tool for this purpose, since the branching ratio 
is essentially dependent on the content of f . For example, if f is a pure ss state, the 
branching ratio is ~ 10 -5 , while if /o is composed of four quarks then the branching ratio 
is expected to be ~ 10~ 4 . 

The strong coupling constants q^k+k- and gf K+K- are among the important hadronic 
parameters entering to the analysis involving and / (980). Indeed, the kaon loop diagrams 
contributing — > / 7 are expected to be in terms of gf K+K~, as well as q^k+k-- The 
coupling constant gf K+K- is studied in light cone QCD sum rules [10] (more about light 
cone QCD sum rules and and its applications can be found in [11, 12]). 

In the present work we calculate the strong coupling constant q^k+k- in light cone 
QCD sum rules method. It should be noted that this constant can be obtained from 
experimental data on meson decays. The goal in the present work is twofold: Firstly, can 
we get new information about the quark content of meson comparing experimental data 
with theoretical results? Secondly, how does light cone QCD work for the asymmetric case, 
i.e., with different Borel mass parameters corresponding to different mass channels? 

The paper is organized as follows. In section 2, we derive sum rules for the q^k+k- 
coupling constant. In section 3, we present our numerical results and conclusion. 



2 Sum rules for g^K+K- coupling constant 

In this section we calculate the strong coupling constant g^x+K- i n light cone QCD sum 
rules. This coupling constant is defined by the following matrix element: 

(K-(g)0°(p,5)|K + (p + g)) , (1) 

where the momentum assignment is specified in brackets and is the polarization vector 
of the meson. In order to calculate the strong coupling constant g^K+K- we consider the 
following correlator function 

iWp,9) = « / d 4 xe^(K( q )\T{4(x)j^(o)}\o) , (2) 

where the quark current = u^^s is the axial vector current and = S7 M s is the 
interpolating current for the meson. 
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The correlator function, in general, can be written in terms of the following five inde- 
pendent invariant functions 

IV G», q) = rii^ + n 2 PnPv + n 3 p^ + n 4 g^ + n 5 q^q u . (3) 

Therefore, our first problem is to choose the kinematical structure. For this aim, we consider 
the phenomenological part of the correlator function. This part can be written as 

(^-(g)0 o (p)|K + (p + g))(^ + (p + g)|J / f|O)(O|^|0(p)) 
V ^ (p 2 ~ m 2 )[(p + q) 2 - m 2 K ] ■ U 

The matrix elements entering Eq. (4) are defined as 

{K + (p + q)\J«\0) = f K (p + q),, 

(0\4\<l>{p)) = rntUe,, . (5) 

Using Eqs. (4) and (5), we get for the physical part 



(p 2 — mV)[(p + q) 2 — m 2 



K. 



(Pn + %) (qu + ^) ■ (6) 



It follows from this expression that the only the structures p^q U) q^q^, q^Pv and p^Pv give 
contribution to the correlator function. In further analysis, we will choose the structure 
p^q u from which the corresponding invariant structure 

jt _ 94>k+k- ficm&U tn\ 
- (P 2 - m 2 )[(p + q) 2 - m 2 K ] ' 1 > 

follows. 

Our next task is the calculation of the correlator function from QCD side. This calcula- 
tion can be carried out by using light cone operator product expansion method, in which we 
work with large momenta, i.e., — p 2 and — (p + q) 2 are both large. The correlator function, 
then, can be calculated as an expansion near to the light cone x 2 ~ 0. The expansion 
involves matrix elements of the nonlocal operators between vacuum and the kaon states, 
i.e., in terms of kaon wave functions with increasing twist. 

After lengthy calculations, we get the following expression for the invariant function 
which is proportional to the structure p^q v 

iv, ( P + 9 f) = ifK /; *.{ ^ + - 4 fli '"»+ G ''"» (i + ^) 

+ ^^r) + ifx £ du u J Vat A [2^.(0,) + <p, + 2<f ± ( ai ) 
+ if K du da 3 dm-^ [2(p ± (oii) - <p\\ + 2(p ± (ai) - <p\\(aii) 

+ » /J [du ( U-l) / o 1 da 3 4 ^ (Q3)fa + , "|| 1+Q3( "- 1)1} 

+ £ du [ da, / o "" da, iF(a ' )bq + + }, (8) 
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where 



and, 



A = m 2 — (p + qu) 2 , 
Ai = m 2 s -[p + q(l + (u-l)a 3 ] 2 , 
A 2 = m 2 s - [p + q(ai + ua 3 )f , 



az rl-t 

dt I dai$>(ai, 1 — ati — t,t) , 
o Jo 



F(a 3 ) = 

F(oii) = - f 1 dt$(t, l-a 3 -t, a 3 ) , 
Jo 

$(ai) = <p\\(ai) + (p±(ai) + <P\\(<Xi) + <p±(ai) ■ 
The functions in Eq. (8) are defined as 

(K(q)\u(x) W A0)\0) = -if K q„ C due iuqx [(p K (u) + x 2 gv 

J 



u 



+ !k (x, - ^\ £ due^ x g 2 (u) , 



{K(q)\u(x)a^ 5 s(0)\0) = i( % x v - Q^)^^ [ due^^u) , 



and 



ru 

G(u) = - g 2 (u)dv . 
Jo 

The matrix elements involving quark-gluon field are determined as 

(K(q)\u(x)^ 5 g s G a p(ux)s(0)\0) = 



I x a q^ \ / x p q^ 

qp [ g a ^ - —^r ] - qp\ gfa 



qx 



qx 



UJb J 



qx 

(K(q)\u(x)^g s G af3 (ux)s(0)\0) 



/ x a q^\ / x p q„ 

qp g a n ] -qpigpn 



qx 



qx 



+ ifK^- (q a xp - qpx a ) f Va t ^(a t )e^ +U ^ , 



(9) 

(10) 

(11) 
(12) 



(13) 
(14) 

(15) 



(16) 



(17) 



where G a p = \e a p pa G pa ', Van = da\da2da 3 5{l — a\ — a 2 — a 3 ). 

The sum rule for g^x+K- is obtained by equating the phenomenological, Eq. (7), and 
theoretical, Eq. (8), parts. 
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In order to suppress the contributions of the continuum and higher states, we perform 
double Borel transformation over the variables — p 2 and — (p + q) 2 on both sides of Eqs. (7) 
and (8), and obtain the following expression for the correlator function 

f K m^ mKK e- m ^l M h- m y M ' = f K e- m V M2 {M 2 VK (u )+4u g 2 (u ) 



- 4[0i(uo) + G 2 (u )} + y^W - 4-^§[<?i(u ) + G 2 (u )} + f ^ 



i / da\ 

u -a 3 



+ 

where 
M 2 



1 /■]-— 0:3 

+ / da 3 / c?q;i 

1— no Juq— Q3 

$(t*i) 2 dF(a 



«3 «3 (itti 




m - ai 2(pj_(ai) + V||(ai) + 2v5±(«i) 



«3 
1 

1-uo 



dec; 



Ffl 



a 3 
a 3 ,0,a 3 ) 



«3 



(18) 



M 2 M\ 
Ml + Ml 



2 ' 



U 



M\ 



M{ + .Uj 



2 ' 



and, 



m 2 K u (l - m ) 



Subtraction of the continuum and higher states is carried out by employing the quark- 
hadron duality, i.e., continuum contribution, which is represented in terms of the spectral 
density obtained from QCD side, by equating it to the one obtained from QCD side, but 
starting from some given threshold. The prescription for subtraction the contribution of 
the continuum in light cone version of the sum rule is proposed in [13] (see also [14]). In [13] 
and in many works, the symmetric point Ml = M| = 2M 2 (i.e., u = 1/2) is considered, 
and then the continuum subtraction is implemented by means of the simple substitution 



-m 2 /M 2 



e -m 2 /M' 2 _ e -s /M 2 



in the leading twist term (in our case leading twist term is the wave function <y?^(-u)). But 
this prescription is not adequate in our case, where the Borel parameters and masses of 
different channels are not equal. In the present work we will follow the analysis given in 
[10], where the prescription for continuum subtraction through use of the Borel parameters 
with different masses in the respective channels is proposed, and properties of the wave 
functions are exploited. Namely, the leading twist-2 wave function can be exploited as a 
power series 



in order to calculate its contribution in the duality region. Here we will neglect the con- 
tinuum subtraction in the higher twist terms altogether, due to their small contribution to 
the theoretical part of the sum rules. Here, we will neglect the continuum subtraction in all 
higher twist terms, due to their small contribution to the theoretical part of the sum rules. 
The final result for the q^kk coupling is given as 



9<t>KK 



m 



e ^l e ^ K /Ml e - m 2 /M 2 i M 2j2 bk ( Ml 

{ k \ M 1 



! _ e -(-o-m2)/M* £ 1 ft 



s - m; 



M 2 



4 



2 \k+l 

s - m z s \ 



+ Au g 2 (u ) - 4[^i(tto) + G 2 (uq)] 



+ 



m 



K 



M\Ml (k + l)\ \ M 2 ) 

Ul 2 ( /■l—uo ruo rl rl—az 

( M o) - 4 T^|[fi , i( M o) + G 2 {uq)\ + / da 3 / da ± + da 3 da x 

M z V JO JO Jl-un Juo-a-i 



- / da 3 / d«i 2 2 — {2(p±(a>i) + ip\\(a>i) + 2ip ± (a>i)) 



+ 



<S>(ai) 2 dF( ai 



a 3 a 3 da\ 



da- 



l—uo 



F'(a 3 ) 

«3 



da. 



F(l - a 3 ,0,a 3 ) 



l—uo 



a 3 



(19) 



where s is the smallest continuum contribution. 



3 Numerical analysis 



In this section we present our numerical calculation on q^kk coupling constant. It follows 
from Eq. (19) that the main input parameters are the kaon wave functions. The theoretical 
framework for their determination is based on an expansion in terms of the matrix elements 
of conformal operators [15]. In particular, for the leading twist -2 wave function ipk(u) 
defined in Eq. (13), the expansion goes into Gegenbauer polynomials: 



Pk(u, aO = 6w(l — u) 

where a 2 (l GeV) = 0.2 [16]. 
Analogously ip a is defined as 



i + EM/^ 2 (2«-i) 

n=l 



(20) 



(p a (u) = 6u(l — U) 



1 



7 2 3 



i - ( -m - ^mw 3 - -p z - -p 2 a 2 ) c 2 3/ \2u - 1) + 



(21) 



where, at the fx = 1 GeV scale r] 3 = 0.015, w 3 = —3, a 2 = 0.2. Here, the factor p = m 2 s /m 2 K 
takes into account the boson mass corrections (see [17]). The twist-4 wave functions ip\\(ai), 
<P±{oti), <P\\((Xi) and <p±(ai), including the meson mass corrections are given as (see [15] and 
[17]) 



h 



h 00 + h 01 a 3 + ^(5o 3 - 3) + • 



(f±(ai) = 30772^0:3(20;! — 1 — a 3 ) 
<P±(cei) = 120m^oi(l — a\ — a 3 )a 3 [ai (l - 2a\ - a 3 ) + 



<P±(<Xi) = -30m x o 3 |/7 00 (l - a 3 ) + h 01 [a 3 (l - a 3 ) + 6«i(l - a x - a 3 ) 
+ h 10 [a 3 (l - a 3 ) - -[a\ + (1 - a x - a 3 f] H 1 , 



<f±(oii) = 120m^oi(l - «i - a 3 )a 3 [voo + ^io(3o3 - 1) + 
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where 



^oo — 


1 


^oo = -^Va , 
3 


hoi = 


7 3 


4 20 


ho = 


7 3 






21 


Vio = 


y^4^4 , 




21 9 


a 10 = 


T^ 4 - 20° 2 



with 774 (// = 1 GeV) = 0.6 and w±{n = 1 GeV) = 0.2 [15, 17]. 

The values of other input parameters appearing in Eq. (19) are: m s = 0.14 GeV [18], 
m K = 0.4937 GeV, = 1.02 GeV. Leptonic decay constant of meson, = 0.234 GeV, 
follows from the experimental result of the — > decay [19]. The threshold So which 

is varied around the value s — 1-1 GeV 2 , is determined from the analysis of two-point 
function sum rules for f K [20]. 

Having all input parameters, we now proceed by carrying out numerical calculation. 
The dependence of g^KK on Borel masses M\ and Mf at two fixed values of s — 1.1 GeV 2 
and s = 1.2 GeV" 2 is presented in Figs. (1) and (2), respectively. According to the QCD 
sum rule method ranges of the auxiliary Borel parameters Mf should be found such that 
the result for q^kk be practically independent of them. 

From these figures we see that, such regions indeed do exist. When M\ and Mf are 
varied in the regions 2 GeV 2 < M\ < 4 GeV 2 and 0.8 GeV 2 < Mf < 1.4 Ge\/ 2 , the result 
for g^xif seems to be independent of the Borel parameters. It should be noted here that, the 
result changes slightly when the continuum threshold is fixed to the value So = 1.2 GeV 2 . 
The final result for g^KK is 

g+KK = 4.9 ± 0.8 . (22) 

At this point, let us discuss sources of the uncertainties. SUf(3) breaking effects in kaon 
distribution amplitudes which we neglected, can play essential role, since we can explore 
wide range of u and hence smoothing the effects of the shape of wave function. Additional 
uncertainty arises from the value of m s . All these factors can cause an uncertainty about 
5-10%. Moreover, the errors coming from the variations in the continuum threshold and 
Borel masses, change the result about 10%. If all these uncertainties are taken into account, 
the resulting error is about 20%, which is quoted in Eq. (22). 

Finally, we would comment that, existing experimental results on — > KK decay 
predicts g$KK — 4.8. So, obviously, we see that our result is quite close to the experimental 
value. Therefore we conclude that the quark content of is ss, and for channels with 
different masses and different Borel parameters, light cone QCD sum rules work quite well. 
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Figure captions 

Fig. (1) The dependence of the coupling constant q^kk on the Borel parameters M\ and 
Mf , at the fixed value s = 1.1 GeV 2 of the continuum threshold. 

Fig. (2) The same as Fig. (1), but at the fixed value s = 1-2 GeV 2 of the continuum 
threshold. 
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Figure 2: 
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